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Abstract 

Using a geometric definition for the lattice Chern-Simons term in even dimen- 
sions, we have studied the distribution of Chern-Simons numbers for the 2d-U(l) 
and the 4d-SU(2) lattice Higgs models. The periodic structure of the distributions is 
preserved in our lattice formulation and has been examined in detail. In both cases 
the finite size effects visible in the distribution of Chern-Simons numbers are well 
accounted for by the Haar measure. Moreover, we find that (Nq S ) grows with the 
spatial volume. We also find numerical evidence that tunneling in 4d is increased 
at high temperature. 



1 Introduction 



Recent interest in the finite temperature Electroweak phase transition has focused 
on the role of tunneling between topologically distinct vacuums. Such tunnelings, 
in a semiclassical approximation, are well known in the continuum and lead to a 
periodic structure in the effective potential During such a tunneling the so- 
called topological charge will change by an integer. In two dimensions there are 
the kinks, which are time independent finite energy solutions of the scalar Higgs 
model or the Sine-Gordon model. The vacuum solutions occur at spatial infinity 
and the tunneling is between <fr(x = — oo) and <p(x = +00), <fi being the Higgs field. 
In the 2d - Abelian gauge Higgs system the tunneling is in Euclidean time and 
governed by time dependent finite energy solutions - the vortices. The tunneling is 
most transparent if one assumes an axial gauge A = 0. The vacuum gauge field at 
temporal infinity becomes pure gauge Ai(t, x) = <f>(t, x)~ 1 di<fi(t, x) and the tunneling 
goes from Ai(t = — 00, x) to Ai(t = +00, x). In four dimensions one has instantons 
in the pure SU(2) gauge theory. This case is very similar to the previous one, upon 
replacing the Higgs field with a proper gauge transformation. With Higgs fields 
included no time dependent finite energy solutions are known. However, assuming 
static fields only, one has a saddle point solution - the sphaleron ||. The sphaleron 
in itself does not provide the tunneling. It is so to speak lying midways between two 
vacuums and we must still imagine a time dependent field configuration interpolating 
among the two vacuums. 

The physical relevance of these tunnelings was pointed out by 't Hooft who found 
that neither baryon nor lepton numbers are conserved in the electroweak theory . 
For the baryon and lepton currents one has 

V? = d,j; = ^trlFM, (1.1) 

where Nf is the number of families of quarks and leptons. While the B — L symme- 
try remains unbroken due to the anomaly cancellation, B + L is no longer conserved. 
This so-called baryon number violation is caused by the nontrivial topological wind- 
ing of the SU(2) gauge fields. The baryon number B changes by an amount 
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In the axial gauge, A = 0, AB is related to the change in the Chern-Simons(CS) 
number Ncs, 

AB = N f [N cs (t 2 ) ~ Ncsih)] (1.3) 

with 

Ncs = - ^ / ^xe^triMdjAk + ^A^)]. (1.4) 

07T J O 

At zero temperature only quantum tunneling via instanton like configurations is 
possible. The rate for such a process is, however, exponentially suppressed. This 
is because the relevant field configurations have an action of the order of the bar- 
rier height 2% /aw, with aw denoting the electroweak coupling constant. At high 
temperatures such an exponential suppression is absent, since tunneling can occur 
classically by thermal fluctuations. Assuming that the temperature is so high that 
only static fields (kinks or sphalerons) are relevant, one can go over to a Hamiltonian 
formulation of the theory and study the evolution of the system in real time via the 
classical Hamiltonian equation of motion or the Langevin equations with a friction 
term. Classical nucleation theory can then be applied to extract information about 
the rate of baryon number violating processes. This has been done for the above 
mentioned lattice models @, |5|, ||, [7j. However, a priori it is not known whether the 
static field approximation is valid and we, therefore have attempted to study the 
full theory in Euclidean time. By varying the time extent of the lattice we can in 
principle control the temperature. We have initiated work in the 2d - Abelian gauge 
Higgs model and the 4d - SU(2) gauge Higgs model ||. 

In two dimensions it is easy to derive an expression for the Chern-Simons num- 
ber, Ncs, on the lattice, it amounts to evaluating the sum of the link angles in a 
Polyakov line. In four dimensions this is much more involved. There are several ge- 
ometric definitions of N C s |], 0' H3> E3- We decided to use the version by Seiberg 



T0|| j which shares the properties of the continuum expression. It is gauge dependent 
and changes by an integer under large gauge transformations. We note that Ncs 
is in general non-integer; only for configurations which behave as pure gauge it is 
indeed an integer. Numerical investigations of topological properties of lattice gauge 
theories are known to be difficult and time consuming, if one attempts to preserve 
the important quantization properties of topological objects also on the lattice. In 
the geometrical approaches an interpolation of the original gauge fields into the ele- 
mentary cells of the d-dimensional lattice is required and a d-dimensional integration 
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over the fields needs to be performed. Starting from an existing program for the 
calculation of topological charges we developed a program for the calculation of 
Chern-Simons numbers, which is fully vectorized and very efficient. 

A central goal in studies of CS-numbers in Higgs models on Euclidean lattices is 
the calculation of the distribution of CS-numbers, P(Ncs), and the determination 
of the temperature dependence of the barrier height between topologically distinct 
vacuums. In the high temperature limit, where static field configurations dominate, 
this is related to the difference between the minima and maxima of the effective 
potential, V(N C s), for the CS-numbers, which can be extracted from the probability 
distribution, 

P(N CS ) = eM-V(Ncs)) 

= J dAd(/>exp(-S(A,<t>)) 

lr. 2 

S(N CS + d 3 xe ijk tr[Ai{dj A k + g4A)D • ( L5 ) 



However, it also is known that even in the continuum at zero temperature, vac- 
uum fluctuations may lead to rather large values of the CS-numbers and may easily 
dominate the above distribution and cover any expected periodic structure of the 
potential expected to arise from topological properties of the field configurations 
H, |^, [j]]. On the lattice one encounters in addition the problem of finding a for- 
mulation for the CS-numbers which correctly reproduces the topological properties 
known in the continuum. It is the purpose of this paper to systematically investigate 
these problems on the lattice. We will study the properties of a lattice version of 
the CS-term suggested by Seiberg WU . 



The outline of the paper is as follows. In Section 2 we discuss topological charge 
and Chern-Simons numbers in the continuum. In Section 3 we give the definitions 
of the lattice topological charge and Chern-Simons term. In Section 4 we derive 
expressions for the distribution of Ncs in two dimensions at strong coupling and 
we present results for the periodic structure. Section 5 is devoted to our results in 
four dimensions, and we describe the algorithm used in our simulations of the SU(2) 
Higgs model. Finally we give our conclusions in Section 6. 



3 



2 Topological charge and the Chern-Simons term 
in the continuum 



We will first consider the 2d - U(l) case. It is well known that one can define a 
gauge invariant topological charge Q: 



Q = ~t~ / dtdxenuFuy £ Z. 

4:71 J M 

= d,A u -d u A^ {2.1) 



The manifold is denoted M and we will assume that its boundary dM is a one 
sphere S 1 . The topological charge density, q, can be written as a divergence of the 
Chern-Simons density K^, 

= ~^vA v . (2.2) 
Under a local gauge transformation, g, the gauge field, A^, changes like 

8A li = \g- 1 d li g , (2.3) 

% 

so that 

SK^-^-e^gg- 1 . (2.4) 
We define the Chern-Simons number Ncs at a given Euclidean time as follows: 



cs 



= [ dxK . (2.5) 

J dM 



We note that iVcs is an integer only for pure gauge (vacuum) configurations and 
is gauge dependent. It changes, however, under a gauge transformation only by an 
integer 



5N, 



cs 



/ dxd 1 gg- 1 eZ. (2.6) 
Ini JdM 



2ixi JdM 

This follows also from homotopy theory using the mapping g : S 1 — > U(l) = S 1 . 
Such mappings are characterized by the homotopy class I1 1 (S' 1 ) G Z. In the tun- 
neling picture vortices will interpolate between two vacuums with integer Ncs — 
n, n + 1. 
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In the 4d - SU(2) case the topological charge Q is 

Q = -ttttt / dtd 3 xe flupa tr[Fp U Fp a ] G Z , 

OZ7T z JM 

Fp V = d^Av — d v An + [An, A v ] . (2.7) 

Here we shall assume that dM = S 3 . Like in two dimensions we can write 

Q 22^2 ^i*vpa^ \Fp V F pi ^\ dpKp 

= -—e^ pa tr[M9 P A a + -A p A a )}. (2.8) 

Under a local gauge transformation, g, the gauge field changes as 

5Ap = g- l \A, + d,}g , (2.9) 

giving 

SK,t= - 2^e^ pa tr[d„g g' 1 d p g g- 1 d a g g' 1 ] 

- ^ep Upa d u tr[d u g g- 1 A a ] . (2.10) 

The (timelike) Chern-Simons number Ncs is then defined exactly as in eq.(2.5), 
and has the same properties. Its gauge variation is an integer (the boundary term 
vanishes) 



5N, 



cs 



^ovp« I <Pxtr[d u gg- l d p gg- l d a gg- l }eZ. (2.11) 

47T ^ JdM 



24tt 2 JdM 

This time the mapping g : S 3 — > SU (2) = S 3 , and the homotopy class is ILj^ 3 ) G Z. 
In this case the instanton will interpolate between two vacuums, assuming the axial 
gauge. 



3 Topological charge and the Chern-Simons term 
on the lattice 

We will now consider the lattice version of the topological charge and the Chern- 



Simons number. We will use a geometric definition given by Seiberg [IIJ. Problems 
with dislocations will be ignored here. It is convenient to start from the definition of 
the topological charge first given by Luscher || . The following considerations apply 
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to 2d or 4d. The manifold M is a torus and we will cover it with cells c(n), where n 
denotes the lattice sites. Let the gauge potential be defined on c(n) and likewise 
A™~^ be defined on c(n — ft). At the faces f(n,fj,) = c(n — p) H c(n), we can relate 
the two potentials by a transition function u„ )/Lt (x), 

A^(x) = *,-»L4"(x) + d u ]v n ^x). (3.1) 

At the corners of the faces the transition function is given by 

Vn^x) = w n ~^x)w n (xy 1 . (3.2) 

Here w n (x) is a parallel transporter, used to gauge fix the links to the complete axial 
gauge in each cell. By interpolation this formula is extended to the whole face. The 
transition function, v n ^(x), defines a bundle, while w n (x) (given on the boundary 
of the cell) is a section of the bundle [ITS] . In Liischers approach M the interpolating 
fields in the interior of the cells are given by 

S n^( x ) = w t {n)S n>ti {x)w t {x)~ 1 , t = n,n- ii . (3.3) 

From now on we shall assume that t = n or t = n — fi and that v ^ \i. In 2d one 
has explicitly 

Sn tli (x) = [w\n)U{n,u)w t {n + v)- 1 \ x , 

w t (x) = [w t (n)w t (n + v)~ 1 ] x w t (n), v ^ /i, t = n,n - fi, 
S n>fi (x) = U(n,v) x . (3.4) 

The approach of Seiberg is quite similar flID| . He does not apply any axial gauge 
fixing, and his function S n>ll depends only on the original fields. However, otherwise 
the interpolation formulas for s* (x) and S njfJi (x) are identical. Although the 4d 
case is more involved the previous statement holds as well. For the Liischer charge 
U one has 

Q L = E^H = E(- 1 ) M (^,m- W.A \q L (n)\< 1/2, 
(-1)^ = —t^l dxsl^xY^sl^x). (3.5) 
Notice that kt tfX is gauge invariant. After a little algebra one finds 

Q L ( n ) = 7T~ l °g[U(n, l)U(n + 1, 2)U{n + 2, l)" 1 ^, 2)" 1 ], \q L {n)\ < 1/2. (3.6) 
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The Seiberg charge jn| is obtained by replacing 

[s^(x),k t J^[S^(x),K t;fl ] , (3.7) 
where K t)tl can now be interpreted as the local Chern-Simons term. 

Q S = 52q S (n) = 52(-inK nt „-K n+ ^), \q s (n)\<l/2, 

n n,fi 

(-lTK ti , = --^ / dxS^{x)~ l d u S t ^x) . (3.8) 

ZTil J f(t,fi) 

Though each K n ^ term is gauge dependent, the charge remains gauge invariant, in 
fact q s (n) = q L (n). This follows from the relation 

q s (n) = q L (n)-q w (n) , 



Q W (n) = —e^ / dxw n (x)- 1 d v w n (x) . (3.9) 

2m Jdc(n) 

The last piece is an integer since we integrate the section over the boundary of the 
cell. By restricting the charge to the interval \q{n)\ < 1/2, it follows that the two 
charge definitions agree. The Chern-Simons term for the 2d-U(l) Higgs model is 
now given by 

^bs = £#«.,,, = ^r£ log tfK,i/) • (3-10) 

The summation is only over the spatial lattice at a fixed Euclidean time, i.e. Ncs = 
Ncs(t)- Under gauge transformations N C s indeed changes by an integer. Moreover, 
it is an integer for pure gauge. 

We can extend these considerations to the four dimensional case. Defining for 
shortness 

K(x) = sl^d^x) , 



one finds 



VUx) = P^M^dvP^M , (3.11) 

Q L = £<?>) = B-im,,- W) , 

n n,fi 

-mt,» = ~^ 2 e, vpa J^^xtr[Sl{x)Sl{x)Sl{x)\ 

+ /, d'xtrlV^SUx)] . (3.12) 



The actual expressions for (x) and Pt+via>{ x ) are gi ven i n re f- 0- In Seibergs 
version we use the same expressions but replace 

K»(x),P?+u,nv(x),k t J [StA x ), p t+i>M x ), K t,n\ ■ (3-13) 

The Seiberg charge will be gauge invariant, if we use the same restriction as in two 
dimensions, and the Chern-Simons term is given by Ncs = J2n s Kn a ,n- I n the naive 
continuum limit a —>■ one finds with U(n,fi) = exp(aA n ^) 



•iyKt^ = -^e^triA^dpA^ + ^A t>p A tt<7 )} . (3.14) 



It is easy to show that 

K t ^ = k t ^ + K^ , (3.15) 

where Kf has the same form as K ttfl , but with S ttfl (x) — > u>*(x). As the section 
transforms like a gauge transformation this implies that the Chern-Simons number 
will change by an integer. For pure gauge fields one finds in addition ^(x) = 1 
and S ttl j,(x) = Pt+v^v(x) at the plaquette. We, therefore, find 

(-l)"# t ,„ = --^p, / d 3 xtr[W:w;w:], W: = w n (x)-%w n (x). (3.16) 

This expression closely resembles the continuum form. After summation over the 
spatial volume Ncs indeed becomes an integer even on finite lattices. 

We finally note that the topological charges, obtained from the two approaches 
discussed above, are related fH| . Using eqs.(3.2) and (3.5) one has for any configu- 
ration 

q s (n) = q L {n)-q w {n) 

qW{n) = ^w/ Sc(fi) ^WW p X] . (3.17) 

Here, q w (n) is the topological charge (integer) of the section and Q L = Q w . Notice, 
that there is no restriction on q L (n) so that q s (n) = q L (n) up to integers. For smooth 
fields like instantons they always agree, while for realistic field configurations this is 
true for almost every cell. 
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4 Chern-Simons number in the strong coupling 
limit, finite size effects 



We will discuss here the distribution of CS-numbers induced by the Haar measure on 
finite lattices, i. e. in the absence of any Boltzmann weight factors in the Euclidean 
path integral over configuration space. We will call this the strong coupling limit. 
For the U(l) Higgs-model in two dimensions it is easy to give an expression for the 
distribution of Chern-Simons numbers in this limit. Using the definition given in 
eq.(3.10) one finds on a lattice with spatial extent n the following recursive relation 
for the density p n (z), of CS-numbers z, 

P ^ Z) = { I ItherwCe 



rmin(^i,z+i) 
i=l"-2 j=l JmaxC- 

Using this recursion formula we find explicitly 



Pn(z) = 11, dz iPl ( Zi )5(J2zj - z) = / dzp n -!(z) . (4.2) 



1 n-l-m / cy \ 

p 2n {z) = ^—^ E (-1)* ^ J ( n - i - \z\f-\ m<\z\<m + 1.(4.3) 
These relations allow us to define moments of the distributions as 

/« 
2 tfc* a p n (z) . (4.4) 
-f 

Some properties of p n (z) and these moments are discussed in the appendix. In 
particular it is obvious that all odd moments vanish and p n (z) / p„(0) — > 1 for n — > oo. 
However, on a lattice of finite extent n, the distribution functions, p n , are close to 
Gaussians, i.e. large values of the CS-numbers are suppressed due to finite size 
effects. This is reflected in the even moments, 

V° = 1 

n 



V A = 

n ' 12 ' 



n n 2 



V 4 = — — + — . (4.5) 
120 48 { 1 
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The leading term for V^ a = (2a — l)!!("P 2 ) a is just coming from a Gaussian 
approximation^. Therefore (Nq S ) grows with the number of links in a timeslice. 
Notice the fourth order cumulant — 3 (7-^) 2 = — 

We will mainly be interested in the gauge invariant, non- integer part of the CS- 
numbers. This can be obtained by summing over the various gauge equivalent sectors 
of the distribution functions. On a lattice with (2n+l) sites in spatial direction these 
are given by (see Appendix eq.(A.l) for even number of lattice sites) 

n 1 

R%n+i(z) = E (z + mr P2n+1 (z + m) , \z\<- . (4.6) 

m=—n 

We show in the appendix that on a lattice of size n the distribution R° n and the 
moments are, in fact, independent of z in the strong coupling limit, and equal 
the moments V%. In particular, this means that R° n {z) = 1 for all z, i.e. the 
distribution of the non-integer part of the CS-numbers is flat in configuration space 
even on finite lattices. 

Let us now discuss some numerical results obtained in the strong coupling regime. 
We used a 2d-U(l) lattice Higgs model. The lattice action is 

s= - -EMM + ^(-«*)-if 

n,fi T n 

with r = 1 for U(l) and r = 2 for SU (2). We studied the distribution of CS-numbers 
for various parameter sets on lattices of size 12 2 and 24 2 . In the left column of Fig.l 
we show three distributions of the CS-numbers obtained from calculations on the 12 2 
lattice at (J3, A) = (2.0, oo) and re = (0.5,2.0,8.0). At each of the re values 100.000 
measurements of the CS-number, defined by eq.(3.10), in each of the 12 time-slices 
has been performed. The distributions for all 12 slices have then been added in 
order to increase the statistics. 

The (3 value has been chosen such that our strong coupling calculations should 
be applicable, whenever the Higgs sector decouples from the gauge field sector. This 
happens at small values of re. Indeed, no periodic structure is visible in the distribu- 
tion for our lowest value of re and the distribution is well described by our analytic 
a The second moment has previously been evaluated in ref. pq| . 
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Table 1: Moments of the distribution of Chern-Simons numbers for the 2d-U(l) 
Higgs model on a 12 2 lattice at (f3, A) = (2.0, oo). The last column shows the strong 
coupling results given in eq.(4.5). 



K 


0.5 


2.0 


8.0 


0.0 


(N 2 CS ) 


0.998(5) 


1.001(8) 


0.970(13) 


1.0 


(N 4 CS ) 


2.886(21) 


2.909(34) 


2.830(90) 


2.9 



result, eq.(4.3). At large values of k, however, the periodic structure becomes trans- 
parent. More and more of the configurations are close to pure gauge, hence the 
clustering of N C s around integers. We note, however, that the peak heights are 
distributed according to the characteristic Gaussian form of the strong coupling 
distribution function. In the right column of Fig.l we have divided the original dis- 
tributions by the strong coupling form of the distribution, pn(z), which reflects the 
non-uniform distributions of CS-numbers in phase space due to finite lattice effects. 
The peak heights are now identical (The large fluctuations at higher Ncs values are 
due to limited statistics in this exponentially suppressed part of the distribution.). 
Although the shape of the distributions changes drastically, the moments of Ncs 
are insensitive to this as can be seen from Table 1. For each k value we find on the 
12 2 lattice (N cs ) ~ 1.00 and (N cs ) ~ 2.90, which is in perfect agreement with the 
exact results Vf 2 = 1 and Vf 2 = 2.9. Likewise for the 24 2 lattice, (N cs ) = 2.09 and 
(N cs ) = 11.80 compared to P| 4 = 2 and V\± = 11.8. We also have checked that 
these results are independent of (3. Hence, for large lattices the data are thus very 
well described by a Gaussian distribution. 

The results presented so far are in accordance with the assumption that the dis- 
tribution of CS-numbers on a lattice of size n a x n T is a product of two probability 
distributions, where one is just given by p„ CT and describes the phase space restric- 
tions for Ncs and the other gives the probability to find a certain non-integer part, 
z, for the CS-number, Ncs = m + z. This part of the distribution is controlled 
by the action and thus contains the relevant physical information. Let us write the 
probability to find CS-number Ncs = m + z as 

Pn a ,n T (N C s) =U natnT (z)p na (m + z) \A <\ » ( 4 - 8 ) 
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Table 2: Moments of the distributions of the non-integer part (z) of Chern-Simons 
numbers, Ncs = m + z, for the 2d-U(l) Higgs model on a 12 2 and 24 2 lattices at 
(P,\) = (2.0,oo). 



K 


(z 2 ) , n = 12 


(z 2 ) , n = 24 


{z 4 ) , n = 12 


(z 4 ) , n = 24 


0.5 


0.0832(1) 


0.0833(1) 


0.0125(1) 


0.0125(1) 


2.0 


0.0662(1) 


0.0805(1) 


0.0091(1) 


0.0119(1) 


8.0 


0.0287(1) 


0.0511(1) 


0.0024(1) 


0.0062(1) 



with p na defined in eq.(4.3) and u> ntrtriT denoting a probability distribution that de- 
pends on the couplings of the Euclidean action as well as the spatial and temporal 
size of the lattice. Using eq.(4.6) and (A. 13) it is then easy to verify that indeed the 
moments, (Ng s ), are independent of (3 and k. 

Let us now discuss the distribution of the non-integer part of the CS-number 
and the corresponding moments. These describe the fluctuations of the CS-number 
around vacuum configurations in a given topological sector. In the continuum theory 
arguments have been given that these vacuum fluctuations are also proportional to 
the spatial volume of the system ||. In Table 2 we give results for the moments of 
the non-integer part of the CS-numbers. We note that for a flat distribution one 
finds (z 2 ) = 1/12 and (z 4 ) = 1/80. These limiting values agree with our results at 
k = 0.5 and are also approached for larger values of k with increasing size of the 
lattice. 

5 Monte-Carlo results in 4d 

Before discussing the results for the 4d-SU(2) Higgs model, it is necessary to give 
some details of the programs developed by us to evaluate K t41 . Because the inter- 
polation in the Seiberg case is done on the original links and not on the gauge fixed 
links, the fields tend to be very rough. The same happens if one wants to evaluate 
the topological charge q w (n) via the section. One of the integrations in eq.(3.12) 
can be done in analytic form, so we are left with a two dimensional integral. These 
integrals must be evaluated carefully, if one wants to extract the periodic structure 
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in the distribution of CS-numbers. On the other hand, if we are only interested in 
the non-integer and therefore gauge invariant part of Ncs (used in the discussion 
of baryon number violation) then it suffices to perform several Landau gauge fixing 
sweeps before doing the integrals. The effect is that Ncs will be shifted by an inte- 
ger, in most cases into the interval — 1/2 < Ncs < 1/2. The advantage is that gauge 
fields on the links become very smooth and the integrals will converge much faster. 
To define a convergence criterion we note that K t jX must change by an integer under 
a gauge transformation. We took various configurations and monitored Ncs, first 
without and then with a Landau gauge fixing sweep. If 5K n ^ was an integer within 
a relative error of e = 10~ 4 , the integrals were accepted. 

We have used the following strategy, which turned out to be very efficient. The 
volume term in eq.(3.12) is the most time consuming part, so we concentrate on 
that. Since the Chern-Simons numbers on each timeslice are hardly correlated we 
can use all of them. We perform a Gauss-Legendre integration with 8x8 points and 
store the results for all the K n ^s in each timeslice. We then repeat this calculation 
with 16 x 16 points and compare the results for each K n ^. If the relative difference 
is less than e the contribution is accepted. Otherwise we collect the K n ^s which 
have not yet converged. In these cells we repeat the integration with 32 x 32 points 
instead, compare with the previous values and eventually repeat the procedure with 
64 x 64 points. This part can be done very efficiently in vectorized form. At this 
point only in a few cells the results for K n ^s have not converged, for these we use a 
library integration routine with interval adaptation. The typical time for evaluation 
of Ncs (without any gauge fixing) in a time slice with 6 3 lattice points is roughly 60 
seconds on the CRAY-YMP. With gauge fixing this is reduced to about 6 seconds. 
The entire program runs with a speed of about (150-200)Mflops. 

As a test of our programs we calculated the time dependence of Ncs on a 4d - 
SU(2) instanton configuration in the complete axial gauge. This is shown in Fig. 2 
for an instanton of core size p = 2, on an 8 3 x 12 lattice. On the first timeslice 
the gauge fields are in the vacuum sector with Ncs — 0. On the last timeslice the 
instanton has mediated the tunneling to the other vacuum sector with Ncs — 1- 
We note that in the middle (t = 6) one finds Ncs = 1/2- The corresponding field 
configuration on this timeslice is the analog to the sphaleron in the Higgs models 
and interpolates between the two vacuums at t — 1 and t = 12. 
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For the SU(2) gauge Higgs model we used the action given in eq.(4.7) with r = 2. 
We have chosen to study the structure of the CS-number distributions for various 
values of k at fixed (/3, A) = (2.25,0.5). For this choice of parameters the phase 
transition between the symmetric and broken phases is known to occur at k ~ 0.27 
IB| . We have performed simulations on a 4 3 x 2 lattice at k = 0.25, 0.30 and 0.40, i. e. 



on both sides of the transition line. For each of the k values we have measured Ncs 
on 2000 gauge field configurations without any gauge fixing. In the left column of 
Fig. 3 we show the CS-number distributions. The periodic structure is obvious even 
for k below k c . With increasing k the CS-numbers on many of the configurations are 
close to an integer, and therefore the gauge field configurations can be interpreted as 
being close to pure gauge. We checked that N C s indeed changed by an integer under 
a gauge transformation, which confirms the quality of our numerical integrations. 
At k = 0.3 we have performed axial as well as Landau gauge fixing. Even after a few 
gauge fixing sweeps nearly all CS-numbers get shifted into the interval [—1/2, 1/2]. 
The distribution after only a single gauge fixing sweep is shown in Fig. 4. With a 
few more Landau gauge fixing sweeps only the peak in the middle would remain. 

For the moments of Ncs we fhid on the 4 3 x 2 lattice at k — 0.30 the values 
(N cs ) = 2.93 and (N cs ) = 24.57. We note that a Gaussian fit with width 2.93 would 
suggest (N cs ) = 24.96. Like in two dimensions the moments are thus very close to 
those of Gaussian distributions. Moreover, we also find for the 4d-SU(2) Higgs model 
that the moments are independent of (/3, k) and increase proportional to the spatial 
volume. We have performed simulations on a 6 3 x 2 lattice (without gauge fixing) at 
(P,X) = (2.25,0.5) and k = (0.30,0.40). For this lattice we find (N cs ) = 9.63 and 
(N cs ) = 274.39. A Gaussian approximation would lead to (N cs ) = 275.26. The 
general structure of the distributions thus is very similar to the 2d-U(l) case and 
we expect that the finite size effects can be eliminated similarly by dividing out the 
distribution induced solely by the Haar measure (strong coupling distribution). We 
approximate this by a Gaussian distribution with the width given by our numerically 
determined value for (N cs ). This is shown in the right column of Fig. 3. The peaks 
are now more or less of equal height. We thus may expect, that similar to the 2d 
case the finite size effects in the distributions drop out in the distribution functions 
uj na ^ nT {z) for the gauge invariant, non-integer part, z, defined in eq.(4.8). We note 
that these can be calculated on gauge fixed configurations, which is computationally 
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much less demanding. 

At finite temperature the tunneling between topologically distinct vacuums should 
become more likely with increasing temperature. This should be reflected in a 
flattening of the distributions, 0J na ,n T { z )i defined in eq.(4.8). In order to see, whether 
this effect is visible in our distributions, we have performed simulations on lattices of 
size 6 3 x n T with n T = 2, 4 and 6. For these lattices we used 50 Landau gauge fixing 
sweeps, to make the integrals converge fast. For both lattice sizes we collected 6000 
Chern-Simons numbers at k — 0.3. The results for uj na ^ nT (z) are shown in Fig.5. 
The left column shows the change in the distributions with varying n T at fixed n a . 
The tendency for a flattening of the distribution at finite temperature is clearly 
visible. We take this as evidence that the system tunnels more often at non-zero 
temperature. 

A major problem for a more quantitative analysis of the distributions at finite 
temperature is caused by the volume dependence of oj nt7t n T (z) itself, which leads to 
a flattening of the distributions with increasing spatial volume. This is seen in the 
right column of Fig.5, where we compare distributions on 4 3 x 2, 6 3 x 2 and 8 3 x 2 
lattices. Clearly the distributions become flatter with increasing spatial lattice size. 
We note, that this effect, caused by the spatial volume, is opposite to the finite 
temperature effect shown in the left column of Fig.5. In that case the distributions 
become more peaked around N C s = 0, although the total 4- volume increases. 

The moments, (z 2 ), of the distributions shown in Fig.5 are summarized in Table 
3. They show that the limiting value for a flat distribution, (z 2 ) = 1/12 is rapidly 
approached with increasing spatial volume as wells as decreasing couplings f3 and/or 
K. As mentioned before the broadening of the distribution with increasing spatial 
volume is not unexpected. It has been shown || that even at zero temperature 
quantum fluctuations of the fields can lead to large values of Ncs- In fact, this 
contribution is proportional to the spatial volume. In the present parameter range 
studied by us such a linear dependence is not visible. We will have to work at 
much large values of the couplings (3 and/or k. This will be necessary in order to 
proceed with a quantitative analysis of finite temperature effects on the CS-number 
distributions. 
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Table 3: The width of the distribution of the non-integer part (z) of Chern-Simons 
numbers, N C s = m + z, for the 4d-SU(2) Higgs model on lattices of size x n T at 
(P, A) = (2.25, 0.5) and various values of k. 





N T 




(z 2 ) 


4 


2 


0.30 


0.028(2) 


6 


2 


0.30 


0.068(2) 


8 


2 


0.30 


0.083(1) 


6 


2 


0.30 


0.068(2) 


6 


4 


0.30 


0.049(2) 


6 


6 


0.30 


0.049(2) 


4 


2 


0.25 


0.057(2) 


4 


2 


0.30 


0.028(2) 


4 


2 


0.40 


0.009(2) 



6 Conclusions 

We have studied the properties of Chern-Simon numbers in two and four dimensional 
Higgs models on Euclidean lattices. We have shown that the definition of the Chern- 
Simons term based on the geometric definition of topological charges preserves the 
basic properties of the continuum expressions. The non-integer part of the Chern- 
Simons number is gauge invariant; the Chern-Simons number changes by an integer 
under gauge transformations. Moreover, we have shown that the effective potential 
for the Chern-Simons number is periodic with maxima of equal height at integer 
values of Ncs, if finite size effects are taken into account properly. 

The distribution of Chern-Simon numbers flattens with increasing temperature. 
This will lead to an increase of baryon number violating processes at high temper- 
ature. At present we could, however, not explore the temperature of the transition 
rates quantitatively as the distributions are still influenced by finite size effects 
caused by quantum fluctuations in the trivial vacuum. These fluctuations will be 
suppressed at larger values of the gauge coupling, which we might have to choose 
for a future quantitative analysis of the Chern-Simons number distributions in the 
vicinity of the symmetry restoring phase transition of the 4d-SU(2) Higgs model 
JT7|] . Another possibility to suppress the contribution from vacuum fluctuation is 
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based on a modification of the cooling method, which also allows to smoothen con- 
figurations and obtain the contribution of classical configurations that extremize 



the Euclidean action [18]. This approach has recently been used in the context of 
the 2d 0(3) nonlinear sigma model |1| and is in an gauge invariant manner also 
applicable to SU(2) gauge theories |20| . It will be interesting to test with our ge- 
ometric formulation whether the Chern-Simons numbers develop a plateau under 
extremization of the Euclidean action. This would allow to extract the distribution 
of Chern-Simons numbers on topologically non-trivial gauge field configurations for 
which the quantum fluctuations are suppressed. 
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Appendix 

We will first show that the following identities are true on lattices with spatial extent 
(2n + 1) and 2n, respectively, 

n y 
E P2n+l(z + m) = 1 , \z\ < -. 
m=—n 
n-1 

E P2n(z + m) = 1, 0<Z<1 . (A.l) 

m=—n 

One has 



n n—m 



(2n)\ E P2n+i(z + m) = E E(-!)M \ {n i - z - m) 2n 



771=1 1 = 

n n—m 



E E(-1)M 2U t l )(n+ 1 -- i + z-mr 



m=l i=0 



+ E(-ir( 2n , + 1 )(^+^-^-N) 2n - (a.2) 



i=0 
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Interchanging the summation over (i, m) in the double sum gives for the first term 



£ El" 1 ) \{n+--i-z-m) 

m=l i=Q \ 1 ) Z 



i=0 m=0 \ / 

The summation over m can now be performed giving for the first term in eq.(A.2) 

£(-ir( 2 f ){n-\-i-z? n ■ (A.3) 
Using now also the recursion relation, 

( 2 7 1 )=( 2 r) + ( i -i) ■ ™ 

in the third term of eq.(A.2), we obtain 



Z " 1 2n \ . 1 



(2n)\ £ p 2n+1 (z + m) = J2(-iy ™ („ - - - i + zf n = (2n)\ . (A.5) 

m=—n i=0 \ / 

This concludes the proof. 

Let us now calculate the moments of this distribution. We consider the integrals 

/n 
dzz a p 2n 
-n 

O n-l n-l-m / n \ , m +i 

- sis_E E (-i)' T / - « - 1* 2 ""' ■ (a.6) 



Interchanging the summation over (i,m) in the double sum and then integrating 
over z yields 



The Stirling numbers of the second kind, 

Sn m = ^±pE(-iy( 7 )<» , (a.8) 

with the recursion 



m! ~ V z 



^E^+n-™-!, 5^ = 1 (A.9) 

1=1 

can be used to write the moments as 

— (2n)!d E(-r-f 2 :!rV^ • (a.io) 
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In the following we shall only consider the cases a = (0, 1, 2, 3, 4). From the recursion 
formula it follows easily 



This leads to 



cm 



s 



rn 



( m + 1 \ 
{ 2 ) > 

3m + 1 / m + 2 



S 



rn+2 - 4 V 3 / 

m(m + 1) / m + 3 



m+3 2 V 4 

m 2 + 5 

48 \ 5 



15m 3 + 30m 2 + 5m - 2 / m + 4 , 



-po 

' 2n 

' 9, 



2n ~~ w ' 



P. 



2 



n 



2n g ' 

Pf„ = , 

PL = + ■ 

We also have checked that the following relations for the "discrete" momenta 
hold as long as the degree of the moment does not exceed the lattice size, 



n ^ 

(z + m) k p 2n+1 (z + m) = Vt+i , k < 2n + 1 , \z\ < -. 

m=~n 
n-1 

(z + m) k p 2n (z + m) = P k n ,k<2n, 0<z<l . (A.13) 
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